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^ ' Abstract 

' The Wheeler-DeWitt equation for the minimally coupled FRW-massive-scalar-field min- 

I isuperspace is written as a two-component Schrodinger equation with an explicitly 'time'- 

, dependent Hamiltonian. This reduces the solution of the Wheeler-DeWitt equation to the 

. eigenvalue problem for a non-relativistic one-dimensional harmonic oscillator and an infinite 

■ series of trivial algebraic equations whose iterative solution is easily found. The solution 

' of these equations yields a mode expansion of the solution of the original Wheeler-DeWitt 

equation. Further analysis of the mode expansion shows that in general the solutions of the 
Wheeler-DeWitt equation for this model are doubly graded, i.e., every solution is a super- 
J> , position of two definite-parity solutions. Moreover, it is shown that the mode expansion of 

[X^ ■ both even and odd-parity solutions is always infinite. It may be terminated artificially to 

' construct approximate solutions. This is demonstrated by working out an explicit example 

which turns out to satisfy DeWitt's boundary condition at initial singularity. 
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1 Introduction 



It is well-known that the Wheeler-DeWitt equation for the minimally coupled FRW-massive- 
scalar-field minisuperspace can be written in the form of the Klein-Gordon equation in (1 + 1) 
dimensions 

^(a,(/)) = 0, (1) 



where a = In a, a is the scale factor of the FRW metric 



ds^ = -dr + a' (t) [dr^ + S^{r) {dO' + sin^ 6 dip' 

sinhr for k = —1 : open universe 
:= { r for k = : flat universe 
sinr for k = \ : closed universe, 

and is a real scalar fleld of mass m. It is also well-known that unlike the massless case exact 
solutions of this equation are not known. The approximate solutions are usually obtained using 
the semiclassical (WKB) methods, 0. 

Recently, I have used a two-component formulation of the Klein-Gordon equation to study the 
geometrical phases induced on a Klein-Gordon fleld due to a dynamical background spacetime 
geometry and electromagnetic or scalar interactions, [|] . In the present article, I shall try to apply 
the same approach for the treatment of the Wheeler-DeWitt equation (|l]). 

2 Two-component Form of the Wheeler-DeWitt Equation 

The idea of the two-component formalism |^ is to transform the Klein-Gordon equation which is 
of second order in time-derivative, into a system of two coupled differential equations which are 
flrst order in the time-derivative. The latter may be viewed as a Schrodinger equation. In this 
way one can apply the well-developed techniques of solving the Schrodinger equation to obtain 
solutions of the Klein-Gordon equation. The application of the two-component formalism for the 
Klein-Gordon equation is described in detail in Ref. Here I only include the results which are 
relevant for the solution of Eq. (|I|). 

"'^The scalar field is assumed to be homogeneous. 

^This equation corresponds to simplest choice of factor ordering. 



d d 
'±_ _^ yfx) 



Eq. ([T|) can be expressed in the form: 

^ ^ = , (2) 

where /i, z/ = 0, 1, {t]^") := diag(— 1, 1), X° := a, X'^ := 0, and 

Hence, a plays the role of 'time.' Furthermore, defining 

"\/2 ^ ■?/; — i-?/^ ^ ' ^ ^ 

one can easily check that Eq. @ is equivalent to the Schrodinger equation 

i-^ = H-^ , (5) 

with the explicitly 'time'-dependent Hamiltonian: 

2 [l-D -l-D I 
where a dot represents differentiation with respect to X^ = a and 

D := V = —— -V . (7) 

The two-component Hamiltonian (Q) has the following remarkable properties: 

— It belongs to the (dynamical) Lie algebra sus{l, 1) ® 1). Here sus{l, 1) is the standard 

representation of sm(1, 1) with the basis 

1 i 

K3 = J3 = -0-3 , Ka = iJa = -<Ja , (o = 1, 2) , 

where are Pauli matrices, and 1) denotes the unitary representation of 1) 

constructed out of the position x and momentum p operators of the one-dimensional quantum 
mechanics, with the basis {(x^ — (x^ + {xp + px)/4}, 0. This is easily seen, 

by writing H in the form H = K_ + K^D, where K± := ± K2, and observing that the 
operator D belongs to sUu{l, 1). 
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It is the square root of the operator in the sense that it satisfies 

D 



D 



D 



(8) 



This is reminiscent of Dirac's treatment of the non-interacting Klein-Gordon equation in 
(3 + l)-dimensional Minkowski space, where the Dirac operator plays the role of the square 
root of the Klein-Gordon operator. 

Next let us examine the eigenvalues En and eigenvectors \E'„ of the Hamiltonian H. It is not 
difficult to show that due to the particular form of H , the eigenvalue equation 



(9) 



leads to the following relations: 



" V2 1-B„ ' • 



= E $ 



(10) 



Therefore the eigenvalue problem for H is reduced to that of D. Substituting Eq. (|^) in the second 
equation in (p!0D, one has: 

r 1 1 _i 



V - UJ 

2 86^ 2 



2 J,2 



11) 



where u := me^° = ma^ and Sn := {E^ + Ke^")/2 = {E^ + Ka^)/2. Eq. (|TTD is the eigenvalue 
equation for a simple harmonic oscillator with frequency uj and unit mass. Hence: 



Ljj(n -) = ma (n + -J 



(12) 
(13) 



where if„ are Hermite polynomials and Nn are normalization constants. The latter are determined 
by the requirement that be normalized with the usual L^{WC) norm. This leads to 

Nn 

Note that and consequently ^n tend to zero as a ^ or a — > oo. 
Substituting Eq. (|T^) in the expression for one has: 





1/4 


m 


1/4 


7r(2"n!)2 




7r(2"n!)2 





E±„ = ±^Juj{2n + 1)- Ke^" = ±y^me3"(2n + 1) - /te"^" = ±a^ma{2n + 1) - ko^ . (14) 
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Thus, for K = — 1 or 0, E^n are real, whereas for k = 1 and a > m, there are eigenvectors 
with imaginary eigenvalue. They correspond to the integers n satisfying m{2n + 1) < a. 

So far, I have given the solution for the eigenvalue problem for the Hamiltonian H. The 
eigenvectors of H, however, do not solve the Schrodinger equation. This is mainly because H 
is explicitly 'time'-dependent. The general solution of the Schrodinger equation for an arbitrary 
explicitly time-dependent Hamiltonian is not known. However, in the cases where the eigenvalue 
problem for the Hamiltonian is solvable, there are techniques based on the adiabatic approxima- 
tion [0, H, ^ which may be used to obtain approximate and occasionally exact solutions. These 
techniques can only be employed if the Hamiltonian is self-adjoint. In order to examine the self- 
adjointness of the Hamiltonian H, one first needs a notion of an inner product on the space Ti. of 
two-component wave functions \E'. 



3 Hilbert Space Structure of Ti 

As a vector space, the space Ti. of two-component wave functions is isomorphic to L^(IR) ©L^(IR). 
Here by the vector space L^(1R)©L^(1R), I mean the vector space inner sum of two copies of L^(IR), 
not the Hilbert space inner sum of these spaces where the sum is also endowed with a canonical 
inner product. Furthermore, in what follows I allow for negative norms and the term Hilbert space 
refers to a complete not necessarily positive-definite inner product space |jlO[. Ti. = L^(IR) ©L^(IR) 



may be viewed as a restriction on the solutions of the Wheeler-DeWitt equation, since for finite 
a the elements of L^(IR) © L^(IR) vanish at 'spatial' infinity, i.e., (p = ^oo. This choice seems, 
however, to be a natural and at the same time a rather general one. It is only restrictive as far 
as the boundary condition at 'spatial' infinity is concerned. It is still quite general, for it does 
not correspond to any fixed boundary conditions at a = ±oo. Furthermore, although Ti has two 
vector subspaces with inner products, an inner product on Ti. is not a priori fixed. 

A Hermitian inner product on 7i may be given by the inner product on L^(IR) and a Hermitian 
two-by-two matrix h. A particularly advantageous choice for h is h = = diag(l, —1), i.e., 

(^, ^) := (^^1^1) - , (15) 

where ( | ) denotes the inner product on L^(IR). This choice, on the one hand, renders H self- 



adjoint. On the other hand, it reduces to the Klein-Gordon inner product^ for the solutions (^) 
of the Schrodinger equation (^. 

Adopting this inner product on 7Y, one can compute 

= {El^ + E^n)Smn . (16) 

Here I have used (|10]) and the fact that are orthonormal. In particular, {'^±n, ^±n) = 
2Re{E±n) = ±2|Re(-E'„)|. Hence, the eigenvectors \['±„ with imaginary eigenvalues (if they exist) 
have zero norm; they are null vectors. This is of course to be expected since the Hamiltonian is 
self-adjoint.0 Furthermore, note that the null eigenvectors are only present for the k = 1 case. 
Another interesting consequence of Eq. ([T6|) is that 

I for En real , , 

I Ezfzn tor En imagmary. 

Hence, if we view Ti. as the linear span of \E'±n, it can be orthogonally decomposed as 

n = n^-^ ®n^'^^ ®n^^\ (is) 

where n^-\ n^^\ and are the linear span of the negative- norm, null, and positive-norm 

eigenvectors. Note that this decomposition is 'time'-dependent (it depends on a = e°), and not 
every null vector belongs to Furthermore, Eq. ([T7|) implies that the basis {^n} (resp. {^_„}) 
of T^^^-* (resp. Ti'-"-*) is orthogonal, whereas TIq can be orthogonally decomposed into the sum of 
two-dimensional subspaces Ti^^ = Span(\E'„ © -n) with n < {a/m — l)/2, namely 

N 

H^^) = ^nf, (19) 

n=l 

where N := [(a/m — l)/2] is the largest integer part of (a/m — l)/2. 

4 Adiabatic Approximation 

Consider the ansatz 

^ = e^^-^ , (20) 



^For more details see Ref. 

^Note that it is the non-positive definiteness of the corresponding norm that allows for imaginary eigenvalues 
and null eigenvectors to exist. 
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for the solution of the Schrodinger equation where Caa is a complex- valued function of a, and 
M = ±n G Plugging pO|) in (||) and making use of the eigenvalue equation and the form 
([T0| ) of the eigenvectors "^j^, one has: 



-(C^ + E^){1 + Ej,) + lEj^] + ^(1 + Ej,)^ 
iCu + Ej,){l - E^) - lE^ $„ + i{l - E^)<i> 



Adding and subtracting both sides of these equations, one finds: 







(21) 



Here the first equation implies that for h ^ Ann '■= ^($n|$n) = 0. This requirement together 
with the second relation in (0) signify the conditions under which the adiabatic approximation 
would yield an exact solution of the Schrodinger equation. Clearly, Ej^ ^ for the case under con- 
sideration. Thus the adiabatic approximation cannot be exact. However, approximate solutions 
may be sought if Ej^ and Ann are small. 

Performing the necessary calculations, one can show that for arbitrary non-negative integers 
h and n: 



4. 



3i 



E 



nin - 1) 4,n-2 - V {n){n - 1) 4,n+2 

3/2 



m a 



2n + 1 - ^ 

m 



where 



AT 



±77, 



(22) 
(23) 



Eq. (^) indicates that the adiabatic approximation is not valid, i.e., the solutions of the form 
(l20| ) do not exist. 

Another consequence of Eq. (^2]) is that in a mode expansion of a solution \1/ of Eq. (1), the 
even modes '^2n- and the odd modes '^2n-+i do not mix in the course of evolution. This suggests 
the existence of a double grading of the set of solutions. In particular one has the even and odd 
solutions: 

oo oo 
^'even = E ^2A'^2Ar , *odd = E ^2aA+1^2a'+1 • (24) 



A'= — CO 
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5 Mode Expansion of the Solutions 



Consider expressing the solution \1/ of the two-component Wheeler-DeWitt equation as a hnear 
combination of the eigenfunctions \l/jv- of the Hamiltonian: 



(25) 



^^= — 00 



Substituting this equation in (|]) and using Eqs. and (|iy) to simphfy the result, one finds 

oo 

^ [(cT„ + ia_nSn)$n + (T„$J = 0, (26) 

n=0 

oo 

[{a_nEn + iOnEl + a_„E,)$, + a.nEn^^ = , (27) 

where 



n=0 



n = 0, 1,2, ■■• . 



Note that in general Z^o 7^ Zq and cr_o 7^ o"o- 

Next take the inner product of both sides of Eqs. ( P^D and ( [27| ) with Making use of Eq. 
and orthonormality of one has: 



&n + iEn.a^n + - ^ (^ + 1 ) (^ + 2) (T„+2 - - \l n{n - 1) a„_2 = 



£;„cr_„ + iE^an + Kcr_„ + - ^/(n + l)(n + 2) En+2(^- 



-(n+2) 
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n(n - 1) K-20--(n-2) 



(28) 
0,(29) 



where o"_(o) = o"_o. 

In order to analyze these equations further, let us solve for a-n from (|28D and substitute the 
result in (^). Doing the same for cr_(„±2)5 one obtains: 



Er, 



(^n + ^\l{n + l){n + 2) (T„+2 - - yn(n - 1) a-„_2 



(30) 



■&n + - -(n' + n + l)]a, + - y'(n + l)(n + 2) (T„+2 



-- ^^(n - 1) an-2 + (-) V ('^ + + 2)(n + 3)(n + 4) ^^+4 + 



(-)yn(n-l)(n-2)(r2-3)a„_4 = 



(31) 



Note that Eq. (|3lD links only cr2ri,'s, with = 0, 1, ■ ■ ■, thus making the analysis of the coefficients 
of non-negative modes independent of those of the negative modes. 



Eqs. (|30D and (|3TD clearly display the splitting of the solutions of the Wheeler-DeWitt equation 
into even and odd solutions. Since these equations are linear and there is no correlation between 
the even and odd solutions, every solution is a linear combination of an even and an odd solution. 
Hence without loss of generality, one can restrict to the analysis of \E'even and \E'odd- 

Before engaging in this analysis, however, I would like to note that the mode expansion (^) 
for the two-component wave function ^, translates into 

oo 

^=Y. ^n^n , (32) 
n=0 

for the one-component wave function. Thus, it is the coefficients of the zero and positive modes 
that determine the solution of the Wheeler-DeWitt equation (|l]). The negative modes are related 
to the derivative of ip with respect to a. Hence, they are also uniquely determined once the zero 
and positive modes are given. 

In the following sections, I shall show that the mode expansion of the solution of the Wheeler- 
DeWitt equation (|ID is always infinite. This is another way of demonstrating the fact that the 
adiabatic approximation, which corresponds to including only a single mode in this expansion, 
always fails. Nevertheless, if one tries to blindly apply the adiabatic approximation or a slightly 
improved version of it in which one sets all a±n to zero, except cr±n, for some n^,, then Eqs. (p^) 
and yield o"„* +E'^an' = 0. This corresponds to the eigenvalue equation for a one-dimensional 
quantum particle of mass 1/2 and potential V{a) := —E"^ = /te^" — {2n + l)me'^". Specifically, 
an* is the zero-energy eigenfunction. For k = 1, V vanishes at a = — oo. It decreases to attain a 
negative minimum as a increases, and then monotonically increases to +oo as a ^ +oo. Hence, 
if one requires the zero-energy eigenfunction to vanish at a = +oo where V = +oo, then the mass 
m must satisfy a peculiar quantization condition. This is analogous to the mass quantization 
advocated in Ref. However, note that the basic hypothesis of this argument, namely the 

validity of the adiabatic approximation or alternatively the presence of only a single mode in the 
mode expansion (|32D, is never satisfied. Hence, the argument leading to the quantization of mass 
does not apply. 
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6 Even Solutions 



Let us first prove: 

Lemma 1: An even solution \I/even = Z^at ^2AA^2Ar ^-5 nontrivial, i.e., \l/even 7^ 0, if and only if 
cr±o + 0. 

In order to prove this statement, one sets n = in Eqs. (12^) and (P^. This leads to 



CTo + i-Eo(T_o + (T2 = , (33) 



.2 , 3^2 



^ocr-o + ^0(^-0 + ^^0t^0 + ^^2(T_2 = 0. (34) 



Therefore, if (T±o = 0, both 02 and cr_2 must vanish. However, by virtue of Eqs. (0) and (|30|) , 
this is sufficient to conclude that cr2j^ = 0, for all a/- G This completes the proof of Lemma 1. 
Furthermore, one can prove 

Lemma 2: Every nontrivial even solution has an infinite mode expansion. 

The even solutions may nevertheless be constructed inductively in terms of a±o. This is done by 
noting that Eqs. (|30| ) and (|3l| ) yield trivial algebraic equations for a±2(n+i) in terms of a±2n- Per- 
forming an induction on n, one can express the coefficients of all the modes as linear combinations 
of a±o and its derivatives. For example, one has 

1 4 

(^2 = --^{-){ao + iEoa_o) , (35) 



3 

o-Q + 2iEo&_o + (2 ~ ^o)'^o + 22^oO--o 



(36) 



It is easily seen that as a result of this construction a2n involves the first n derivatives of (Tq, first 
n — 1 derivatives of cr_o, and the first n — 1 — j derivatives of E2j, for j = 0, 1, ■ ■ ■ , n — 1. 

Since according to Eq. (|23|), En diverges at the zero of En, demanding regularity of would 
impose certain conditions on cr±Q. In fact, in view of the appearance of E^ in Eq. ([3T|), which is 
analytic for every a, one can easily show that the only source of possible singularity is the term 
proportional to Eq in the expression for (T2. Thus, in order to ensure the regularity of a2n, one 
must require that at the zero of Eo, i.e., a = ln(m/K), ct-qEq and all its derivatives must vanish. 
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In other words there must be an open neighborhood O of a = ln(m/K), in which cr_o(a;) is of the 
form Eo{a)f{a) for some smooth (C°°) function /. In particular a^Q{ln{m/ n)) must vanish. Note 
that for the flat (k = 0) and open (k = —1) cases, this condition is void and a2n are regular if and 
only if a±o are regular. 

Finally note that termination of the mode expansion for a nontrivial even solution is equivalent 
to vanishing of at least three consecutive cr2nS. This would impose three functional conditions on 
(T-to which are in general incompatible. This proves Lemma 2. 

7 Odd Solutions 

The structure of the odd solutions is similar to that of the even solutions. The role of a±Q is played 
by cT-t-i. In particular, one has 

Lemma 3: An odd solution "^odd = Saa ^2aa+i^2aa+i is nontrivial, if and only if <J±i : = 
Zi ± Z_i ^ 0. 



This is a direct consequence of Eqs. ( |28D and ( pi]) which also provide an inductive construction of 
the mode coefficients o'±(^2n+i) in terms of a±i. For example one has: 

= -^i^){&i+tE^a_^) , (37) 
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(3^ 



A natural consequence of this inductive construction is 

Lemma 4: Every nontrivial odd solution has an infinite mode expansion. 

The regularity condition on the odd solutions is analogous to the even case. One can again 
show that the only source of possible singularity is the term Eicj-i in the expression for as. Hence 
the regularity of (J2n+i and therefore ip is ensured by requiring that in the vicinity of the zero of 
El, namely a = ln(3m//t), (T_i must be of the form Ei{a)g{a), for some smooth function g. 

The occurrence of two free functions of a in the construction of even and odd solutions is in 
complete agreement with the general structure of the Wheeler-DeWitt /Klein- Gordon equation (|^) 
whose solutions depend on two boundary or two initial conditions. The latter are given by two 
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functions of a which in the above procedure can be used to determine a±o or a±i depending on 
whether the solution is even or odd. In general the solution will be the sum of an even and an 
odd part. In this case the boundary (initial) conditions would determine both a±o and a±i. Once 
these functions are determined the above iterative construction yields the mode expansion of the 
general solution. 



8 Real Solutions 

Probably the best-known proposal for the solution of the Wheeler-DeWitt equation is that of 



Hartle and Hawking ||T2| who suggested the path integral expression ip oc J exp(— /) for the ground 
state wave function. Here / stands for the Euclidean action. Hartle-Hawking wave function is by 
construction real. 

It is not difficult to see that the reality of the one-component solution of the Wheeler-DeWitt 



equation implies the reality of the coefficient functions cr„ in its mode expansion (|3^). This is 
because are already real. In view of this observation and Eq. (|28|), En<y-n must be imaginary. 
However, as discussed in section 2, for the open and flat universes (k = —1,0), is always real 
so that cr_„ must be imaginary. For the close universe (k = 1) En is real for a < m{2n + 1) and 
imaginary otherwise. Hence in this case (/« = 1) the opposite holds for cr_„, i.e., it is imaginary 
for a < m{2n + 1) and real for a > m{2n + 1). 

These conditions can also be stated as the reality condition on (Tq, cti, (T2, and which also 
generate the solutions. Therefore, they are rather trivial. 



9 Approximate Solutions 

As shown in the preceding sections, the mode expansion of the solutions never terminates. There- 
fore, one cannot use it to construct an exact solution unless one can actually flnd the generic 
expression for the mode coefficients and sum the series. This is an almost impossible task. How- 
ever, one can terminate the mode expansion artificially and obtain an infinite class of approximate 
solutions, ip ~ ipn, := Y^n=o '^n^n- The domain of reliability of this solutions will in general depend 
on the boundary (initial) conditions, alternatively the functions a±o and a±i. 
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In this section, I shall examine the simple choice (T_o = (^±1 = and ctq = 1. Clearly, this 
corresponds to a real even solution whose mode coefficients can be obtained using Eq. (|3T|) . For 
example for the closed FRW case where k = 1, one has 

(72 = 0, a4 = (^)^/'(^ + a^-ma='), = (^)(3m - 4a) , 



4 

(J8 -- 



243 



W^(4a^ - 40ma^ + 36m^a^ + 324a^ - 288ma^ + 189) 



16 

aio = ^(96a^ - lOOOma^ + 7A4m'^a^ + 3604a - 2655m)a^ . 

10935^7^ ^ 



Substituting these relations in the expression for i/jn and making use of the identity |]T3 



[n/2] 

H„{x) = 2"x" + J2 (-l)'2"~^(2z - 1) 



n 
2i 



to calculate the mode functions one can compute ipo = ip2, V^4, V'e, V's and V'lo- Figures 1 to 5 
illustrate three-dimensional plots of these functions. Figures 6 to 9 are plots of the restriction of 
ipQ, ips and tpio to the lines a = 0.1, 1, 10, 100. These figures show that for small a the finite mode 
approximation is quite reliable. In fact, as depicted in Fig. 6, for a = 0.1, the graphs of ips and 
ipio coincide. For larger a one may improve the approximation by including the contribution of 
higher modes. 

Note also that initially V'n are non-zero for a wide range of values of 0. However, as a increases 
the width of this region, which is symmetric with respect to the line = 0, shrinks. This is easily 
seen from Figures 6 to 9. Note also that for large a, the solutions picks about = and the 
approximate solutions become less reliable. As a ^ oo, all the approximate solutions seem to 
vanish except at = 0. 



10 Discussion and Conclusion 

Probably the most important aspect of the two-component formulation of the Wheeler-DeWitt 
equation is the fact that the construction of the solutions in this formulation does not, a priori, 
require fixing a particular inner product or boundary conditions. The freedom in the choice of 
the inner product is refiected as the independence of the mode expansion from the choice of the 
Hermitian matrix h which together with the inner product on L^(IR) determines the inner 
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product on H = L^(]R) © L^(]R). As far as the boundary conditions are concerned, besides the 
imphcit restriction of vanishing boundary condition at 'spatial' infinity, the analysis presented 
here applies generally. 

An interesting consequence of the application of the two- component formalism is a specific 
double grading of the solutions of the Wheeler-DeWitt equation. This allows one to write down 
the most general solution of this equation as a sum of definite-parity (even and odd) solutions 
whose structure can be easily understood. The two definite-parity parts of the solution do not 
mix during the evolution. In particular, if the initial wave function has definite parity, its parity 
remains intact throughout the evolution. 

Unfortunately, the mode expansion of the solutions never terminates. Thus, a closed expres- 
sion for an exact solution cannot be found in this way. However, approximate solutions can be 
constructed by artificially terminating the mode expansion. The domain of applicability of such 
approximations depend generally on the boundary (initial) conditions. 

An alternative approximation is the adiabatic approximation which is, however, not valid. The 
ansatz of the adiabatic approximation, i.e., ip ~ cT„(a)$„(a, 0), is reminiscent of the application 
of the Born-Oppenheimer approximation in quantum cosmology originally advocated by Banks 



1^ . In this approach one views the gravitational degrees of freedom as the 'slow' variables 
and the matter degrees of freedom as the 'fast' variables and considers solutions of the form 
"0 = V'gravity(a)'?/'matter(a, (^)- The rcsults of this article indicate that for ■^/'matter = $n such an 
approximation cannot be valid. 

Another notable implication of the results of this article is regarding the implicit boundary 
condition at the initial singularity a = 0. If one assumes that the coefficient functions (T„ are 
regular at a = 0, then as a consequence of the appearance of a^^^ in the expression for the 
one-component wave-function ip vanishes at a = 0. This is much in the spirit of an argument of 
DeWitt [T^ . Here it follows as a direct consequence of a regularity condition (which may or may 
not hold). As seen from Figures 1-5, the choice cr_o = <^±i = and ctq = 1 leads to approximate 
solutions which are regular at a = 0. They clearly satisfy DeWitt 's boundary conditions. 

A problem which is not addressed in this article is that of the convergence of the mode expan- 
sion in Ti. This is a difficult problem and I have not been able to find the necessary and sufficient 
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conditions under which the mode expansion converges. In fact, a rigorous convergence analysis 
requires a choice of a positive-definite inner product (norm) on Ti. However, as demonstrated in 
the preceding section, there are regions of the minisuperspace in which the contribution from the 
higher order modes is smalL This may be viewed as an indication of possible pointwise asymp- 
totic convergence of the mode expansion in this regions. This argument is however far from being 
satisfactory and further investigation of this point is necessary. 
Finally, I wish to emphasize the following points: 

— The fact that the two-component Hamiltonian belongs to the dynamical Lie algebra 
SMs(l, 1) ® su„(l, 1) suggests a further investigation of the solution of the corresponding 
Schrodinger equation using the group theoretical methods such as those employed in solving 
the Schrodinger equation for time-dependent harmonic oscillators, |jl6|. Recall that the 



Hamiltonian for a time-dependent harmonic oscillator belongs to S'Uu(l, 1), P]. 

The two-component approach to Wheeler-DeWitt equation has a wider domain of applicabil- 
ity than the simple minisuperspace model studied in this article. For example, the addition 
of an interaction potential for the scalar field would only complicate the defining equation 
for and the inductive construction of the mode expansion. The general structure of the 
method would however remain intact. In particular, the 0-dependence of the solutions would 
still be given by $„. More generally, for any other minisuperspace, in the gauge in which 



the shift vector corresponding to DeWitt's supermetric []T5| vanishes, the two-component 
Hamiltonian is of the form (|^), where D is some second order elliptic operator. Therefore, 
even in this case the developments reported in this article are essentially applicable. 



Note: After the completion of this work, Refs. IT^ 113] were brought to my attention where 
the authors use two different two-component formulations of the Wheeler-DeWitt equation. The 
methods described in these articles differ substantially from the present article's. I have also 
noticed a recent paper by Embacher [|l^] who defines the operator D of Eq. (|^) as the (mat- 
ter) energy operator and obtains the mode expansion (|32D in the context of the one-component 
Wheeler-DeWitt equation. He does not, however, obtain the equations satisfied by the mode 
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coefficients. Instead, he focuses on the semiclassical and adiabatic approximations. The two- 
component approach developed in this paper justifies the choice of the mode functions made by 
Embacher in a very natural way. It also makes it possible to obtain the structure of the mode 
expansion by iteratively solving for the mode coefficients. 
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4 

Figure 1: Plot of 'i/'o = "^2 = J2n=o (^n^n, for a±i = cr_o = 0, o"o = 1 and m = 10 
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Figure 6: Plot of the restriction of %pQ (dotted line), ip^ (dashed line) and (full line) to the 
a = 0.1 line. Note that the graphs of ip^ and ipiQ coincide for this value of a. This shows the 
reliability of the finite mode approximation. Note also the range of the values of (j) over which -0^ 
are non-zero. 
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Figure 7: Plot of the restriction of ipQ (dotted line), ip^ (dashed line) and V'lo (full line) to the 
a — 1 line. Note the range of the values of (f) over which -0^ are non-zero. 
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Figure 8: Plot of the restriction of ips (dotted line), and ■0io (full line) to the a — 10 line. Note 
the range of the values of over which ipn are non-zero. 




Figure 9: Plot of the restriction of ips (dotted line), and ■0io (full hne) to the a — 100 line. Note 
the range of the values of over which ipn are non-zero. 
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